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Abstract 

Nucleoli resonance information is often obtained from fitting hadro dynamical 
model calculations to data, where model parameters such as resonance coupling 
constants are the free parameters in the fitting procedure. For reactions with a 
limited data set, such as p{'^,K~^)K, complications in the extraction of reliable 
N* information occur not only through theoretical uncertainties, but also due to 
technical difficulties in the fitting procedure. In this article we outline a fitting 
strategy based on a genetic algorithm and illustrate the kind of ambiguities which 
can arise. 
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Understanding the composite, quark- 
gluon structure of hadrons is a major 
topic in intermediate energy nuclear 
physics. One way to investigate this 
system is to study its excitation spec- 
trum, as it is a direct consequence 
of the hadron's composite character. 
At present, the experimental identi- 
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fication of the full excitation spec- 
trum is far from complete, and most 
of our knowledge is based on analy- 
ses of pion induced or pion produc- 
tion reactions. 

A complete understanding of the 
nucleon resonances can only be 
achieved by a full, coupled-channel 
analysis which incorporates reactions 
involving all possible initial and final 
states to which nucleon resonances 
can couple [1,2,3]. Meson production 
reactions such as — vrA^, irnN, 
r]N, uN, KA, i^S, . . ., should be in- 
cluded, as well as Compton scatter- 



Preprint submitted to Elsevier Science 



8 February 2008 



ing and meson-induced production 
reactions. It is clear that a simultane- 
ous description of all such reactions 
is an enormous task. 

When the (photoinduced) meson 
production processes are treated at 
a hadronic level, the resonances are 
described by effective fields with in- 
dependent quantities such as mass, 
strong decay widths and photo he- 
licity amplitudes. Usually one deter- 
mines these resonance parameters 
by fitting model calculations to the 
available data. The fitting procedures 
are therefore a necessary step in the 
extraction of resonance information 
from experimental measurements. 
However, by taking into account all 
known resonances in a full coupled- 
channel calculation, one would re- 
quire a few hundred parameters to be 
determined. As a first step, a single- 
channel analysis of one particular re- 
action is essential to identify general 
features, and to simplify the problem 
by reducing it to one dependent on 
only a few tens of parameters. 

The search for an optimum fit to 
the data, by a model depending on 
what is still a large number of pa- 
rameters, is not straightforward. The 
situation is considerably worsened 
when reaction channels such as uN, 
KA or KH are examined, since only 
limited data sets are available. In 
these cases, one may encounter the 
situation where model parameters 
are not sufficiently constrained by 
the restricted data sets and simi- 
lar qualities-of-fit are obtained with 
different sets of parameter values. 
In this paper, we highlight the is- 
sue by presenting a fitting procedure 



based on a genetic algorithm (GA) 
and commenting on some of the re- 
lated ambiguities which arise in the 
complex fitting processes. We have 
chosen the ^(7, K'^)A reaction as an 
example with a limited amount of 
published data, and where we can 
restrict the number of free parame- 
ters by employing a single-channel 
model. 

In our model, the ^(7, K^)A reac- 
tion dynamics is described in terms 
of hadronic degrees of freedom by 
adopting an effective Lagrangian ap- 
proach [4]. The tree level Feynman 
diagrams contain the usual Born 
terms, the A'*(892) and Ki(1270) t- 
channel mesons and two hyperon res- 
onances 5oi(1800) and Poi(1810) in 
the w-channel. All these terms consti- 
tute the so-called background. In the 
s-channel, the nucleon resonances 
5n(1650), Fii(1710), Fi3(1720) and 
Di3(1895) are included. Note that 
the Di3(1895) is a candidate for a 
"missing" nucleon resonance [5] . The 
finite extension of the meson-baryon 
vertices is implemented by the use 
of hadronic form factors [6,7]. In 
Ref. [4] we stressed the difficulties 
associated with parameterizing the 
background diagrams in ^(7, K~^)A 
calculations and presented results for 
three plausible background schemes. 
Subsequent work [8] showed that the 
p{e, e'K^)A process is highly selec- 
tive with respect to viable choices 
for dealing with the background di- 
agrams. The model used here is the 
only one which we found to repro- 
duce simultaneously the ^(7, K^)A 
and p{e, e'K^)A data. 

All the extracted resonance parame- 
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ters G[N*] product of an elec- 

tromagnetic and a strong coupling. 
A description of the various types of 
resonance parameters G'[A^*], their 
normalization and their connection 
to the Lagrangian structure is given 
in Ref . [4] . To minimize the number of 
free parameters, one overall hadronic 
form factor cutoff mass for all Born 
terms and one for all resonance di- 
agrams are introduced. Those cou- 
pling constants and cutoff masses 
constitute the free parameters of the 
model. In the present hadrodynam- 
ical model, this amounts to 22 free 
parameters which have to be deter- 
mined by the data. We have used the 
SAPHIR data set [9] which contain 
24 total cross section, 90 differential 
cross section and 12 recoil polariza- 
tion asymmetry points, in the pho- 
ton energy range from threshold up 
to 2 GeV. The fits are performed by 
minimizing x^- Finding the global 
minimum in a 22-dimensional 
parameter space is not trivial since 
although some of the parameters are 
constrained to a specific range by 
physics arguments, in principle they 
are all unknown quantities. 

In order to minimize the function 
and extract A^* information from 
the ]?(7, K~^)A data, we have used 
a genetic algorithm. In this strat- 
egy, a set of solutions is randomly 
generated. Each solution p is an 
encoding of trial values of the free 
parameters, and is used to evaluate 
the (p) function which determines 
its "fitness" . The population is then 
"evolved" in a manner analogous to 
biological evolution, either by mix- 
ing values ( "crossover" ) or by chang- 
ing values at random ("mutation"). 
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Fig. 1. Model calculations for the en- 
ergy distribution of the total p{-y, K^)A 
cross section. The curves correspond to 
the 38 sohitions with < 3. The data 
are from Ref. [9]. 

The population gradually migrates 
to one or more better points in pa- 
rameter space. Whilst these optima 
are not guaranteed to be global, GA 
research [10,11] has shown that "rea- 
sonable" solutions can be found very 
efficiently. 

The output of the GA is the indi- 
vidual with the best fitness after a 
prescribed number of fitness function 
evaluations. The parameters associ- 
ated with this individual are then 
used as input to the Minuit [12] pack- 
age which employs a variation of the 
Davidon-Fletcher-Powell variable- 
metric minimization algorithm. In 
this way we take advantage of the 
wide search capabilities of the GA, 
but use the strength of a traditional 
optimizer to obtain final solutions. 

To test whether global optima were 
being found, 50 GA plus Minuit 
calculations (run concurrently on a 
Linux-based PC farm) were carried 
out. Of these calculations, 38 re- 
sulted in Minuit converging with a 
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Fig. 2. Model calculations as in Fig. 1 
but for the angular distribution of 
the differential p(7, K~^)K cross section. 
The data are from Ref. [9]. 

below 3.00. This was thought to 
be a reasonable value, since our pre- 
vious work [4] had shown values 
of just under 3.00 to be achievable. 
The range of the fitness of these so- 
lutions^ was 2.45 <x^< 2.93, indi- 
cating a reasonable goodness-of-fit to 
the available data for each solution. 
This is apparent in Figs. 1, 2, and 3 
where the model calculations for the 
total and differential cross section, 
and the recoil polarization asym- 
metry are plotted against SAPHIR 
data. 

Given this similarity in the descrip- 



Prom now on, we refer to "solutions" 
as meaning those in the converged set 
of 38. 
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Fig. 3. Model calculations as in Fig. 1 
but for the angular distribution of the 
p{'^,K~^)K recoil polarization asymme- 
try. The data are from Ref. [9]. 

tion of the available data, it is inter- 
esting to investigate how close these 
solutions of the minimization proce- 
dure are to each other in parameter 
space. In Fig. 4, resonance parame- 
ters for all the solutions are plotted. 
At a first glance, it appears that the 
resulting surface in the free pa- 
rameter space is extremely erratic 
and has a large number of local min- 
ima. The solutions for G[Pii(1710)] 
show that the value could be any- 
thing within a wide range. On the 
other hand, the 5'ii(1650) coupling 
constant exhibits clustering in one re- 
gion, indicating that the value must 
be within a reasonably small range. 
Of great interest is the clustering ex- 
hibited in the G^^^ [P13 (1720)] cou- 
pling constant. Here, there are three 
distinct clusters plus one isolated 
point. Furthermore, two clusters ap- 
pear in the -D13 parameters. We are 
therefore unable to pin down clear- 
cut values for all A^* parameters. 
This clearly illustrates the shortcom- 
ings in extracting reliable resonance 
parameters from the available lim- 
ited p(7, K^)K data. 

To study how the coupling con- 
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Fig. 4. Values of the coupling constants of the four N* resonances 5ii(1650), 
Pii(1710), Pi3(1720) and L»i3(1895) in the reaction dynamics of p{-f,K+)A. They 
correspond to the best solutions of the minimization procedure. The four clusters 
are labeled by A T, * and #. 



stants are correlated, we have as- 
signed each of the four clusters in the 
G(^) [Pi3 (1720)] parameter a distinct 
symbol. In Fig. 4, there appears to 
be a significant correlation between 
the Pi3 parameter and the D13 cou- 
plings. In the -D13, the A and ▼ clus- 
ters form one subset, whilst the * 
and + form another. Since the role of 
the Di3(1895) in this p{-f,K+)A re- 
action is still under debate, we have 
split the solutions into two subsets 
according to this clustering in the 
parameters. The characteristics 
of the two subsets are summarized in 
Table 1. 

As a further study, we show model 
predictions for the two subsets for 
some unmeasured observables. In 
Figs. 5 and 6 this is done for the pho- 
ton beam asymmetry S, the beam- 
recoil asymmetry and the beam- 



Subset 


1 


2 


ix') 


2.56 


2.69 




0.13 


0.12 


best 


2.45 


2.50 


number of points 


9 


29 


cluster labels 


▲ T 


★ # 



Table 1 



Characteristics of the two subsets local- 
ized in the solutions of the minimization 
procedure. The averaged (x^), its stan- 
dard deviation and the best x^ are given 
together with the number of points in 
the set and the corresponding labels in 
Fig 4. 

target asymmetry E. Since none of 
these asymmetries for the ^(7, K~^)A 
reaction are available yet, they do 
not represent constraints in the min- 
imization procedure. Within each of 
the two subsets, the predicted kaon- 
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angle dependence of E, Ox and E 
can be seen to follow similar trends 
(the dashed exception in Fig. 5 is the 
solitary T solution). 

More important discrepancies are 
observed when comparing the asym- 
metry distributions of the two sub- 
sets. The main difference appears 
to be a result of the value of the 
Di3(1895) resonance parameters. It 
is interesting to note that the value 
of the G [Pu (1710)] parameter ap- 
pears to have only a small influence 
in the asymmetries. Reliable extrac- 
tion of Pii (1710) information from 
p{j, K~^)A data is therefore likely to 
be difficult, even when more (dou- 
ble) polarization asymmetries are 
available. The model predictions in 
Fig. 6 show that also the differences 
in the P13 parameters between the * 
and # solutions do not translate to 
significant differences in the plotted 
asymmetries. On the other hand, the 
calculation with the "odd" behavior 
in Fig. 5 is the solitary ▼ solution. 
Since this solution deviates from 
the A cluster most strongly in the 
P13 parameters, this shows a non- 
negligible influence of the Pi3(1720) 
resonance on the observables. 

The calculations show that the 
present data can not discriminate 
between two main regions in param- 
eters space. The clear differences in 
the asymmetries between the two 
subsets will require only modest 
statistical accuracy to distinguish 
experimentally. Therefore, the up- 
coming data for the photon beam 
asymmetry [13,14] and the double 
polarized asymmetries [15] will be 
vital in constraining the values of 
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Fig. 5. Model predictions with the so- 
lutions of subset 1 (A and ▼) for the 

angular distribution at various photon 
energies of the ^(7, K~^)A photon beam 
asymmetry (S), the beam-recoil asym- 
metry (Ox) and the beam-target sym- 
metry (E). 
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Fig. 6. Model calculations as in Fig. 5 
but with the subset 2 (* and #). 
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resonance coupling constants and 
understanding the reaction mecha- 
nism underlying K^K production. 

In conclusion, we have developed a 
fitting strategy based on a genetic al- 
gorithm in order to extract resonance 
information from ^(7, K^)K data. 
The advantages of the genetic algo- 
rithm are that the available parame- 
ter space is extensively sampled, and 
that the large number of calculations 
allows one to investigate any ambi- 
guities in the solutions. The success 
of this strategy perhaps points to its 
use in problems with similar or larger 
numbers of free parameters. 

The solutions obtained from the 
minimization runs result in model 
calculations in close agreement with 
the available data set. However, from 
a detailed exploration of the param- 
eter values, it appears that not all 
solutions are contained in the same 
region of parameter space. For some 
of the variables, large deviations ap- 
pear and no well defined value can 
be pinned down. This is particu- 
larly true for the Pu (1710) coupling 
constant. For other A^* parameters, 
clustering in a specific area is ob- 
served. One stable cluster appears 
for the 5*11 (1650) parameter and for 
the Pi3 (1720) and D13 (1895), a few 
isolated groups of solutions are iden- 
tified. On the basis of these clusters 
in the A^* parameters, we were able 
to select two different subsets of so- 
lutions of the minimization. So- 
lutions within each subset produce 
more or less identical predictions 
for unmeasured polarization observ- 
ables, but large differences were ob- 
served when comparing the two sets. 



The calculations indicate that the 
extraction of reliable resonance infor- 
mation from a limited data set is by 
no means trivial, and has to be han- 
dled with the greatest of care. An im- 
portant conclusion is that when us- 
ing theoretical models to analyze the 
excitation spectrum of the nucleon 
based on a limited data set, the ac- 
curacy of predictions will be greatly 
affected by uncertainties stemming 
from the fitting process. Hence, em- 
ploying a strategy which enables a 
thorough exploration of parameter 
space is essential. We have shown 
that the upcoming measurements of 
specific polarization asymmetries are 
essential to reduce the uncertainties 
resulting from the present limited 
data set. 
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